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INFRARED  TRANSMISSION  THROUGH  SCREENING  SMOKES: 
EXPERIMENTAL  CONSIDERATIONS 


I.  INTRODUCTION. 

In  order  to  evaluate  the  effect  of  tactical  screening  smokes  on  infrared  transmission, 
it  is  necessary  to  understand  the  complex  interactions  among  many  variables.  The  transmission 
depends  upon  the  bulk  properties  of  the  material  (e.g.,  index  of  refraction),  as  well  as  the 
particle-size  distribution,  concentration,  and  pathlength.  The  measured  or  apparent  transmission 
can  be  quite  different  from  true  transmission  because  several  simple  basic  facts  are  often 
overlooked.  Therefore,  it  becomes  quite  difficult  to  compare  data  from  different  laboratories  or 
from  field  tests  because  the  experimental  methodology  is  different  at  each  location  although  data 
from  any  one  particular  laboratory  may  appear  self-consistent.  The  variability  in  the  grouped 
data  (all  laboratories)  may  be  extremely  high  due  to  the  methodology.  In  fact,  the  data  in  each 
laboratory  may  be  so  repeatable  that  each  claims  the  other  is  in  error. 

In  principle,  it  is  possible  to  calculate  the  transmission  of  the  smoke  if  the 
particle-size  distribution,  concentration,  and  pathlength,  and  the  complex  index  of  refraction  are 
known.  But  these  parameters  are  not  always  known  precisely,  and  one  resorts  to  experimentation 
to  define  them.  The  experiment  becomes  that  of  introducing  a  smoke  with  unknown  infrared 
properties  between  the  target  and  the  detector.  The  ratio  of  the  signal  received  with  smoke  to 
that  without  smoke  is  taken  as  the  transmission. 

In  addition  to  transmissometers,  the  effects  of  screening  smokes  on  thermal  imaging 
systems  is  of  interest.  The  experiment  is  similar  to  the  transmissometer  test.  A  smoke  is 
introduced  until  the  target  can  no  longer  be  perceived  by  an  observer.  The  amount  of  smoke 
required  is  a  measure  of  the  smoke’s  obscuring  power.  A  second  type  experiment  consists  of 
introducing  a  “standard”  smoke  between  several  different  imaging  systems  and  a  target.  The 
concentration  is  increased  until  the  target  has  disappeared  on  several  imaging  systems.  At  this 
point,  certain  conclusions  are  drawn  about  the  smoke’s  screening  effectiveness,  but  the  results 
may  be  a  measure  of  the  system’s  behavior  and  the  observer’s  ability  rather  than  the  smoke’s 
effectiveness. 

Presented  in  sections  II  through  VII  are  the  basic  principles  that  the  experimentalist 
should  consider  before  designing  the  experiment.  Although  Beer’s  transmission  law  is  appropriate 
for  most  situations,  scattered  light  can  increase  the  detected  signal. 

As  shown  in  section  II,  by  decreasing  the  detector’s  field  of  view  (FOV),  the 
scattered  light  can  be  minimized  so  that  the  measured  transmission  will  be  representative  of  the 
true  transmission.  When  using  Beer’s  law,  it  is  tacitly  assumed  that  the  particle-size  distribution  is 
not  changing.  Any  dynamic  process  such  as  coagulation  or  sedimentation  can  significantly  affect 
the  applicability  of  Beer’s  law. 

For  broadband  detectors,  the  measured  transmission  depends  upon  the  atmospheric 
transmission,  target  temperature,  detector  spectral  response,  and  the  smoke’s  spectral 
transmission.  It  is  impossible  to  separate  these  factors  from  transmission  data  unless  every 
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parameter  is  carefully  measured  and  specified  (section  III).  In  the  infrared,  the  cloud  can 
attenuate  the  target  signature,  emit  radiation  if  hot,  and  reflect  radiation  from  other  sources.  In 
section  IV,  a  detailed  analysis  is  given  on  how  to  separate  these  various  detected  signals. 

Thermal  imaging  systems  operate  differently  from  transmissometers.  For  a  target  to 
disappear  from  the  screen,  its  signature  must  be  below  the  equivalent  noise  level.  This  noise  level 
depends  upon  the  target  angular  subtense  (section  V).  Furthermore,  there  are  two  distinctly 
different  types  of  thermal  imaging  systems.  How  each  of  these  operates  and  how  a  smoke  may 
obscure  a  target  for  each  is  given  in  section  VI.  In  addition  to  attenuating  the  target  signature,  a 
smoke  may  contain  random  hot  spots  which  lead  to  thermal  clutter.  Clutter  can  confuse  the 
observer  and  lead  to  false  target  detection.  Since  random  clutter  is  difficult  to  analyze,  a 
simplified  case  (large  hot  cloud  obscuring  a  portion  of  the  FOV)  is  presented  in  section  VII. 

Finally,  all  of  the  possible  errors  that  may  occur  during  experiments  are  shown  in 
three  flow  charts. 

II,  LASER  TRANSMISSION  AND  SCATTERING. 

Using  Beer’s  law  (also  called  Bouguer’s  and  Beer-Lambert  law),  it  is  possible  to 
calculate  the  transmission  of  an  aerosol  when  the  mass  extinction  coefficient  a,  the 
concentration  C,  and  the  pathlength  L  are  known: 

T=eXP['^LaC(£)d2]  0) 

Although  C  and  L  are  physical  quantities  that  can  be  easily  changed,  a  is  a  function  of  the 
material  and  the  particle-size  distribution  and  shape.  Assuming  spherical  particles,  a  can  be 
calculated  from  the  Mie  scattering  theory: 

3  /  ,  r  \  N(r) 

a  =  4fi  7rQ(m'l)Tdr  ® 

where  Q  is  the  Mie  scattering  efficiency  factor,  m  is  the  complex  index  of  refraction,  r  is  the 
radius  of  aerosol  particle,  X  is  the  wavelength,  and  N(r)  is  the  number  density  size  distribution. 
The  reason  that  a  is  usually  obtained  experimentally  is  that  often  neither  m  nor  N(r)  is  known. 
a  can  be  calculated  using  equation  1  with  the  experimentally  obtained  values  T,  C,  and  L.  From 
equation  2,  it  follows  that,  if  the  size  distribution  changes,  then  a  will  change  also.  Any  dynamic 
process,  such  as  coagulation,  sedimentation,  or  introduction  of  a  new  aerosol  with  a  different 
size  distribution,  will  change  N(r).  For  hygroscopic  smokes,  such  as  phosphorus  and  HC,  the  size 
distribution  will  depend  upon  the  relative  humidity  and,  hence,  a  will  be  a  function  of  relative 
humidity. 


In  recent  experiments,  fog  oil  smoke  was  introduced  into  a  large  chamber  (1.2  by 
3.6  by  1.5  m  high)  by  a  vaporization  condensation  method.  A  laser  transmissometer  operating  at 
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X  =  0.6328  jum  measured  the  transmission  (L  =  1.2  m)  as  a  function  of  time,  as  shown  in  figure  1. 
The  concentration  was  simultaneously  obtained,  as  shown  in  figure  2.  This  behavior  is  rather 
unusual  because  Beer’s  law  states  that,  if  the  concentration  is  constant,  as  shown  in  the  first 
20  minutes  of  the  experiment,  then  the  transmission  should  also  be  constant. 


Figure  1 .  Transmission  of  a  Fog  Oil  Smoke  as  a  Function  of  Time 

Data  were  obtained  with  a  HeNe  laser  transmissometer 
(X  =  0.6328  /im).  Chamber  size  was  1 .2  by  3.6  by 
1.5  m  high.  Laser  pathlength  was  1.2  m.  Since  the 
transmission  is  changing  with  time,  some  type  of 
dynamic  process  is  occurring. 
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Figure  2.  Mass  Concentration  of  a  Fog  Oil  Smoke  as  a  Function  of  Time 

Data  were  obtained  by  gravimetric  analysis  of  glass  fiber  filters.  The  aerosol  was  aspirated  for  1  minute 
for  each  sample.  The  aerosol  is  stable  for  the  first  20  minutes  before  sedimentation  takes  place. 


Fortunately,  the  answer  lies  in  the  size  distribution.  The  distribution  was  obtained 
with  cascade  impactors  at  time  t  =  0.0  min  and  t  =  51  min.  As  shown  in  figure  3,  the  total 
number  of  particles  greater  than  2.0  jum  actually  increased,  indicating  that  significant  coagulation 
has  occurred. 


Figure  3.  Particle-Size  Distribution  Obtained  at 
t  =  0  and  t  =  51  Minutes  for  a  Fog  Oil  Smoke 

Data  were  obtained  from  cascade  impactors.  This  shift 
in  the  mode  diameter  is  characteristic  of  coagulation. 
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Using  the  data  in  figures  1  and  2,  the  extinction  coefficient  was  calculated  and  is 
plotted  in  figure  4.  We  see  that  a  is  a  function  of  time  because  N(r)  is  a  function  of  time  also. 


Figure  4.  Calculated  Extinction  Coefficient  as  a  Function 
of  Time  for  a  Fog  Oil  Smoke 

The  transmission  data  shown  in  figure  1 ,  the 
mass  concentration  data  shown  in  figure  2, 
and  the  pathlength  L  =  1.2  m  were  used  in  Beer’s 
law.  Since  a  is  a  function  of  time,  the  age  of 
the  cloud  is  significant  when  specifying  a. 


As  a  general  rule,  coagulation  becomes  significant  when  the  number  density  exceeds 
10f  particles/cm  ^  and  the  diameter  is  less  than  1pm.  Sedimentation,  on  the  other  hand, 
dominates  when  D  >  1  pm  in  a  chamber. 

The  extinction  coefficient  derived  from  Mie  theory  calculations  assumes  that  any 
light  scattered  out  of  the  incident  beam  contributes  to  total  extinction.  These  calculations 
generally  assume  that  the  incident  beam  is  a  plane  parallel  wave  of  infinite  extent  and  that  the 
detector  FOV  is  essentially  zero.  Beer’s  law  is  always  valid  for  zero  FOV  detectors.  In  reality, 
the  detector  FOV  is  finite  and  the  total  amount  of  light  entering  the  detector  will  depend  upon 
how  much  scattering  has  taken  place.  Middleton  *  has  shown  that  even  for  the  single  scattering 
case  the  amount  of  detected  scattered  light  increases  significantly  as  the  detector  FOV  increases. 
His  results  (figure  5)  are  not  directly  applicable  to  laser  transmissometers  but  they  do  point  out 
the  importance  of  exactly  knowing  the  limitations  of  the  experimental  setup.  Zuev2  did  not  give 
experimental  details,  but  he  showed  that  Beer’s  law  is  valid  down  to  T=10'13  (figure  6). 
Although  our  fog  oil  smoke  experiments  show  significant  coagulation  at  r  =  6,  Zuev  does  not 
mention  any  dynamic  changes,  even  up  to  t  =  30. 
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For  dense  aerosols,  multiple  scattering  can  occur.  Second-order  forward  scattering 
has  been  heated  theoretically.^  The  results  indicate  that,  for  typical  polydisperse  naturally 
occurring  aerosols  (Deirmendjian  haze  M,  H,  and  C3),  a  correction  factor  can  be  applied  to  the 
measured  transmission  to  obtain  the  true  extinction.  The  factor 

ameasured  an 

R  '  ~^T =  =F  (3) 

appears  in  Beer’s  law  as 

Tmeas  =  e_a^CL  =  e'RaTCL  (4) 

As  shown  in  figure  7,  this  factor  depends  upon  the  receiver  FOV  and  the  size  distribution.  This 
correction  factor  has  been  formulated  in  a  slightly  different  form4  to  account  for  all 
narrow-angle  forward  scattered  light: 


>j>  —  g~(uCL  "t  D) 


(5) 


where  D  is  a  function  of  the  size  parameter  (r/X),  receiver  FOV,  and  the  divergence  of  the  light 
beam.  The  values  R  and  D  have  not  been  experimentally  verified.  Mooradian  et  al.$  have  shown 
that  the  detected  signal  is  a  function  of  receiver  FOV  for  naturally  occurring  fogs.  (See  figure  8.) 

The  effects  of  multiple  scattering  are  not  limited  to  transmission  measurements.  The 
effects  of  second-order  scattering  on  laser-induced  detection  and  ranging  (LIDAR)  returns  were 
considered  theoretically  by  Eloranta.^  It  has  been  shown  experimentally  ^  that  the  backscatter 
signal  increases  as  the  receiver  FOV  increases.  (See  figure  9.) 

Thus,  the  measured  extinction  depends  upon  the  particle-size  distribution, 
concentration,  and  pathlength,  and  the  receiver  FOV.  Unless  the  precise  relationship  among  these 
parameters  is  known,  it  is  quite  likely  that  the  measured  a  will  be  different  from  the  true  a. 
Therefore,  laboratories  with  different  FOV’s  may  obtain  different  extinction  coefficients  for  the 
same  smoke  generated  under  identical  environmental  conditions.  Any  dynamic  process  that  can 
alter  the  size  distribution  will  affect  a.  The  value  of  a  obtained  will  depend  on  the  age  of  the 
aerosol.  Furthermore,  since  it  may  take  several  minutes  to  generate  a  smoke  and  to  uniformly 
disperse  it  within  a  chamber,  the  time  defined  as  t  =  0  is  rather  arbitrary. 

III.  BROADBAND  DETECTORS:  SPECTRAL  RESPONSE  CONSIDERATION. 

The  extinction  coefficient,  as  shown  in  equation  1,  is,  in  general,  a  function  of 
wavelength.  For  broadband  detectors,  the  measured  extinction  is  a  function  of  the  source 
temperature  and  the  detector  spectral  response,  as  well  as  the  spectral  transmission  of  the  smoke. 
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RECEIVER  FIELD  OF  VIEW  (MRAD) 

Figure  9.  Detected  Backscatter  as  a  Function  of  Receiver  Field  of  View  (from  Reference  7) 


For  dense  artificial  water  droplet  clouds  as  aC  increases,  the  amount  of  backscatter 
light  incfeases.  An  argon  laser  (X  =  0.5 145  jum)  was  used.  In  the  limiting  case  of 
very  dense  aerosols,  the  amount  of  backscatter  will  be  independent  of  optical 
depth.  (See  curve  2  in  figure  6.) 
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The  total  flux  detected  depends  upon  the  geometry  of  the  scenario.  Let  us  assume 
that  this  is  fixed  and  that  the  detector’s  geometry  is  included  in  the  calibration  of  the 
instrument.  The  flux  reaching  a  detector  from  a  target  with  emissivity  ey  at  a  single  wavelength 
is  given  by 


FOT  “  raeTRT  +  La 

where  ra  is  the  transmission  of  the  intervening  atmosphere,  Ry  is  the  blackbody  power  emitted 
by  the  target,  and  La  is  the  radiance  of  the  atmosphere.  When  an  emissive  smoke  is  introduced, 
the  flux  reaching  the  detector  becomes 

fST  =  WTRT  +  La  +  ra  Rs  (7) 

where  rg  is  the  transmission  of  the  smoke,  La'  is  the  radiance  modified  by  the  smoke,  and  ra'  is 
the  transmission  of  the  atmosphere  between  the  smoke  and  detector.  Rs  is  a  combination  of  the 
flux  emitted  by  the  smoke  and  the  flux  reflected  off  the  cloud  from  external  sources.  This  is  the 
emission-reflection  term. 

The  current  signal  generated  by  these  fluxes  is  given  by 

I  =  SF  (8) 

where  S  is  the  sensitivity  of  the  detector.  For  broadband  detectors,  the  detector  integrates  over 
the  wavelength  of  interest  so  that 


I  =  /  SFdX  (9) 

JX 

All  the  parameters  La,  R,  e,  and  r  are  in  general  a  function  of  wavelength.  With  synchronous 
detection,  the  source  is  modulated  at  some  frequency  f  and  the  detector  electronics  are  sensitive 
to  this  frequency  only.  The  path  radiance  and  smoke  emission  are  dc  components  and  therefore 
are  not  measured.  We  assume  that  the  smoke  emission  does  not  saturate  the  detector  and  that  it 
is  uniform  over  the  entire  FOV.  The  measured  ratio  of  smoke  to  no-smoke  currents  becomes 


<T>  = 


‘ST 

!OT 


f SrsraeTRTdX 
f STa6yRydX 


(10) 


where  <r>  is  the  mathematical  average  with  respect  to  the  function  SraeyRy.  It  is  the  spectral 
averaged  transmission  and  depends  upon  the  spectral  character  of  S,  ra,  Ry,  and  r. 
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The  effect  of  spectral  mismatch  between  the  smoke  and  the  detector  only  is  shown 
in  the  following  exaggerated  example.  Consider  four  different  hypothetical  systems,  all  of  which 
are  nominally  classified  as  8-  to  14-pm  radiometer-transmissometer  systems.  The  spectral 
sensitivity  of  each  is  shown  in  the  table  and  depicted  in  figure  10.  Assume  that  the  geometry  is 
identical  for  each  detector  and  that  scattering  is  not  a  factor,  as  described  in  section  II. 

Table.  Spectral  Sensitivity  of  Four  “8-  to  14-pm”  Broadband  Radiometer  Systems 


Detector 

Spectral  response* 

pm 

1 

8-14 

2 

11-14 

3 

8-12.5 

4 

8-11 

*  Relative  response  is  1 .0  in  these  regions  and  zero  elsewhere. 
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Figure  10.  Spectral  Response  of  Four  Hypothetical  Detectors 


The  relationship  between  these  detectors  and  a 
hypothetical  smoke  is  shown. 
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Consider  a  screening  agent  whose  transmission  is  80%  from  8  to  1 1  /im  and  1 00% 
from  11  to  14  jum  for  unit  pathlength  and  unit  concentration.  As  the  smoke  concentration 
increases,  the  transmission  from  11  to  14jum  stays  at  100%,  and  the  transmission  from  8  to 
1 1  fxm  approaches  zero.  Thus,  as  far  as  detector  2  is  concerned,  the  smoke  is  transparent.  Since 
detector  1  integrates  over  the  entire  region,  the  transmission  approaches  50%  and,  with 
detector  4,  the  transmission  goes  to  zero.  Detector  3  shows  some  intermediate  value.  In  figure  11, 
the  expected  transmission  for  each  of  these  detectors  is  plotted  as  a  function  of  concentration. 
The  values  were  calculated  with  equation  1 0  and  by  letting  R j  =  ej  =  ra  =  1  • 


Figure  1 1.  Calculated  Transmission  as  a  Function  of  Concentration. 


The  extremely  large  difference  in  transmission  between 
detector  2  and  detector  4  illustrates  why  spectral 
responses  must  be  considered  before  collecting  data. 

A  second  way  of  comparing  these  systems  is  to  calculate  the  apparent  extinction  as 
a  function  of  concentration  (figure  12).  This  representation  is  usually  the  method  by  which 
smokes  are  characterized  experimentally.  Suppose  then  a  laboratory  had  the  8-  to  14-jum 
transmissometer  system  with  detector  1 .  After  performing  a  series  of  experiments  at  several 
different  concentrations  (say  C-3,  4,  5,  and  6)  and  performing  a  least-square  fit  to  the 
experimental  data,  the  laboratory  would  obtain  an  extinction  coefficient  of  about  0.05.  Another 
laboratory  with  another  “8-  to  14-jum”  system  with  detector  2  would  say  the  material  was 
transparent. 
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Figure  12.  Calculated  Extinction  Coefficient  as  a  Function  of  Concentration 


The  extinction  coefficient,  calculated  from  transmission  data,  for 
the  four  hypothetical  detectors  can  vary  dramatically  depending 
upon  the  spectral  response  of  the  detector. 

Although  this  example  was  exaggerated,  the  shape  of  the  hypothetical  smoke 
transmission  is  somewhat  similar  to  phosphorus  smoke.8  It  was  pointed  out  by  the  author 
(appendix  A)  that  different  spectral  response  detectors  can  lead  to  different  measured 
transmissions  for  the  same  concentration  and  pathlength. 

The  effect  is  not  limited  to  the  spectral  mismatch  of  the  smoke  and  detector.  The 
spectral  emission  of  the  target  (source)  will  also  affect  the  measured  transmission  or  extinction 
coefficient.  Assuming  that  source  is  a  blackbody,  we  have  calculated  the  expected  transmission 
for  phosphoric  acid  smoke  as  a  function  of  source  temperature  for  a  typical  HgCdTe  8-  to 
14-/im  detector  (figure  13).  We  have  included  the  effects  of  the  atmosphere  (relative  humidity, 
72%;  pathlength,  400  m). 

Thus  it  is  easy  to  see  that  there  is  no  simple  way  of  obtaining  the  smoke 
transmission  with  a  broadband  detector  because  the  measured  transmission  depends  upon  the 
target  temperature,  spectral  response  of  the  detector,  and  the  atmospheric  transmission.  The 
atmospheric  transmission  is  a  function  of  pathlength,  relative  humidity,  temperature,  and  local 
naturally  occurring  aerosols. 

IV.  AEROSOL  EXTINCTION  EMISSION-REFLECTION  MEASUREMENTS. 

In  section  III,  it  was  assumed  that  the  signal  reaching  the  detector  was  modified 
only  by  the  aerosol  extinction.  In  this  section,  we  consider  three  possibilities:  (l)the  smoke 
may  be  hotter  than  ambient  and,  therefore,  may  emit  radiation;  (2)  the  cloud  may  reflect 


20 


- 1 - 1 - 

'NOISSIIAISNVBl  3MOIAIS 


O) 


CO 


o 
-  o 
CO 


1 


X 

(D  o 

ts  D 

<D 

Oh  q 
X  w 

pq 


CO 


3 

.£P 

E 


ro 


21 


Atmospheric  pathlength  =  400  m,  RH  -  72%,  and  CL  - 


radiation  from  external  sources;  and  (3)  the  atmosphere  itself  may  produce  radiation.  Note  that, 
in  field  tests,  the  sun  and  earth  are  the  external  sources;  whereas,  in  the  laboratory,  the  chamber 
walls  are  the  external  sources.  With  care,  it  is  possible  to  eliminate  these  effects. 

There  are  two  general  methods  available  for  separating  extinction  from  emission, 
reflection,  and  path  radiance.  The  first  is  synchronous  detection,  and  the  second  is  the  use  of 
two  detectors;  one  that  looks  at  the  target  and  the  other  that  looks  at  the  background.  The 
second  method  is  equivalent  to  having  a  single  detector  that  alternately  looks  at  the  target  and 
background  and  is  representative  of  all  thermal  imaging  systems. 

Synchronous  detection  was  discussed  in  section  III,  and  equation  10  gives  the 
measured  transmission. 

If  two  light  choppers  are  used,  so  that  the  source  is  modulated  at  frequency  fj  and 
the  total  flux  reaching  the  detector  is  modulated  at  f2,  then  we  can  extract  both  the 
transmission  and  reflection-emission  properties  of  the  smoke.  The  detector  requires  two  tuned 
circuits:  one  sensitive  to  fj  and  another  sensitive  to  f2.  Using  equation  10,  we  obtain  <r>  from 
the  circuit  sensitive  to  fj.  The  ratio  of  no  smoke  to  smoke  for  the  circuit  sensitive  to  f2  is 


lST 
JOT  , 

h 


J  S^s7  aeTRT  +  V  +  ra,Rs^dX 

J'si.T  aeTRT  + 


(11) 


It  is  possible  to  separate  the  terms  if  several  reasonable  assumptions  are  made.  If  the 
atmosphere  is  at  the  same  temperature  as  the  background,  then  the  path  radiance  can  be 
neglected  and 


!ST 

*OT 


=  <T> 


/sv 


Rsdx 


fi  +  ~1 


OT 


(12) 


If  narrow  band  detectors  are  used,  then  the  emission-reflection  term  is  obtained 
directly.  Therefore,  if  we  had  a  scanning  spectrophotometer,  we  could  obtain  the  entire  spectral 
emission-reflection  from  the  cloud. 

If  synchronous  detection  is  not  used,  it  is  easy  to  see  from  equations  6  and  7  that 
there  is  no  way  in  general  of  separating  emission  from  transmission  by  a  single  measurement  with 
broadband  detectors.  These  effects  can  be  minimized  by  judiciously  choosing  wavelengths, 
matching  the  receiver  FOV  and  aperture  to  the  source  size  and  beam  divergence,  and  using 
selective  bandwidth  filters.  For  laser  transmissometers,  extremely  narrow  bandwidth  filters  can  be 
used;  and,  if  external  sources  are  minimized,  transmission  measurements  can  yield  good  data 
provided  that  all  the  precautions  listed  in  section  II  are  considered. 
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Let  us  consider  the  case  where  a  second  detector  is  allowed  to  look  at  the 
background  (or  alternately  one  detector  scans  the  target  and  background  as  with  a  thermal 
imaging  system).  The  flux  reaching  the  detector  due  to  the  background  is 

FOB  =  raeBRB  +  La  (13) 

and  with  an  emissive-reflective  smoke  this  becomes 


FSB  “  TsTaeBRB  +  La'  +  ra,Rs  (14) 

If  we  measure  the  difference  between  the  target  and  background,  the  ratio  of  smoke  to 
no-smoke  signals  becomes  (assuming  La  =  La ) 


<T>  = 


JST  ~  JSB 

!OT  ”  !OB 


y*STsTa(eTRT  _  eBRB^dX 

^^S(£'j’R'j'  —  CgRg)dX 


(15) 


Thus,  for  a  two-detector  system  sensitive  only  to  flux  differentials,  the  detected 
signal  does  not  depend  upon  the  path  radiance  or  the  emissive-reflective  properties  of  the  smoke. 
The  primary  assumption  is  that  the  smoke  is  homogeneous  and  covers  both  the  target  and 
background  uniformly.  For  smokes  whose  transmission  is  wavelength  independent,  then 
ts  =  <rf>.  Note  that  the  average  transmission  measured  with  a  two-detector  system  is  not  the 
same  as  the  average  measured  with  synchronous  detection.  Both  methods  yield  the  true  smoke 
transmission  only  if  narrow  band  detectors  are  used. 

As  shown  in  all  of  the  above  equations,  at  a  certain  wavelength  the  current 
generated  by  the  detector  is  proportional  to  the  flux  reaching  the  detector  (equation  8).  The 
total  current  generated  is  the  integrated  effect  of  the  spectral  response  of  the  detector  and  the 
spectral  emission  of  the  source  (equation  9).  Since  the  source  is  usually  a  blackbody,  the 
detector  response  is  calibrated  against  the  blackbody  temperature.  Note  that  the  detector  is  not 
sensitive  to  the  temperature  directly  but  only  to  the  flux  reaching  it.  The  relationship  between 
photon  emission  and  temperature  is  given  by  Planck’s  blackbody  equation.  Since  this  is  not  a 
linear  relationship,  the  calibration  of  a  detector  is  not  linear.  (See  figure  14.)  For  convenience, 
the  output  is  usually  linearized  about  ambient  temperature  so  that 

1  =  I0  (1  +  K  AT)  (16) 

where  IQ  is  the  current  generated  by  a  blackbody  source  operating  at  ambient  temperature.  The 
difficulty  with  this  relationship  is  that  the  current  is  related  to  the  flux  by  equation  9. 
Introduction  of  a  smoke  with  spectral  character  will  alter  the  calibration  (the  constant  K). 
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Figure  14.  Typical  Calibration  Curve  for  an  AGA  Thermovision  System 

The  output  is  calibrated  as  a  function  of  the  blackbody 
temperature. 

As  shown  in  appendix  B,  a  more  fruitful  analysis  which  indicates  the  magnitude  of 
errors  introduced  is  to  linearize  the  target  output  so  that 

R(Tq  +  AT)  =  R(T0)  [1  +aAT]  (17) 

where  R(TQ)  is  the  ambient  blackbody  flux  evaluated  at  temperature  TQ  and  at  a  single 
wavelength  midpoint  in  the  region  of  interest.  Using  this  relationship  in  equations  10  and  15  and 
if  ej  =  eg,  it  can  easily  be  shown  that 


<T>  =  <T>  = 


(18) 


Although  this  result  is  what  would  be  expected  if  equation  16  is  used,  the  important  difference 
is  that  mathematically  we  have  linearized  the  flux  reaching  the  detector  rather  than  the 
calibration  of  the  system.  Smoke  will  attenuate  only  flux  not  the  calibration.  The  AT  analysis 
conforms  to  convention  but  can  introduce  significant  errors.  The  magnitude  of  the  errors  is 
similar  to  that  shown  in  section  III.  The  use  of  the  expansion  also  introduces  errors  but  it 
conceptually  makes  sense. 
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The  effects  of  atmospheric  transmission  is  readily  apparent  in  all  of  the  equations 
shown.  Since  pathlengths  are  generally  very  short  in  the  laboratory,  ra  is  often  neglected. 
However,  in  the  field,  the  measured  transmission  is  affected  by  ra  and  cannot  be  neglected.  This 
again  emphasizes  the  fact  that  broadband  transmission  measurements  are  confounded  by  the 
atmosphere,  the  target  temperature,  and  the  detector  spectral  sensitivity.  Even  though  emission 
and  reflection  can  be  eliminated,  broadband  detectors  at  best  can  yield  the  spectral  averaged 
transmission  <r>. 

V.  IMAGING  SYSTEMS:  MINIMUM  RESOLVABLE  TEMPERATURE  CONSIDERATIONS. 

In  section  III,  it  was  assumed  that  the  broadband  detector  was  calibrated  so  that  the 
output  could  be  displayed  on  some  type  of  voltmeter.  The  imaging  systems,  in  addition  to 
having  the  same  limitations  as  discussed  in  the  preceding  sections,  also  rely  upon  the  observer’s 
evaluation  of  the  image.  As  mentioned  when  evaluating  the  effectiveness  of  a  smoke,  the  smoke 
concentration  is  increased  until  an  observer  can  no  longer  detect  the  target.  At  this  point,  the 
target  signature  is  below  the  internal  noise  of  the  entire  detector-electronics-human  observer 
system. 


The  noise  level  is  referred  back  to  the  input  as  an  equivalent  temperature  and  is 
specified  as  the  minimum  resolvable  temperature  (MRT)  above  ambient  for  a  “standard”  observer 
and  is  expressed  as  a  temperature  differential  AT  above  ambient.  The  MRT  is  a  function  of  the 
angular  subtense  a  of  the  target  and  is  plotted  for  two  systems  in  figure  15.  These  systems  are 
identical  in  the  sense  that  for  all  targets  larger  than  ac  they  exhibit  the  same  MRT.  At  another 
target  angular  subtense  aj,  the  MRT’s  of  the  two  systems  are  given  by  ATj^jj  and  ATj^ 
respectively.  Assume  that  a  target  is  presented  before  these  two  systems  and  it  is  AT  above  the 
ambient.  If  an  absorbing  smoke  is  placed  between  the  detectors  and  the  target,  then,  since  AT  j 
is  larger  than  AT2,  more  smoke  is  needed  to  bring  the  target  signature  below  the  MRT  for 
system  1  than  for  system  2. 


TARGET  ANGULAR  SUBTENSE 

Figure  15.  Minimum  Resolvable  Temperature  of  Two  Hypothetical  Imaging 
Systems  as  a  Function  of  Target  Angular  Subtense 
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Consider  instead  that,  during  two  different  tests,  only  system  1  is  used  and  the 
target  sizes  are  slightly  different  (aj  and  02).  Since  AT j  is  greater  than  AT3,  more  smoke  is 
needed  to  obscure  target  size  aj.  During  field  smoke  trials,  actual  military  targets  (tanks,  trucks, 
etc.)  are  used.  The  equivalent  AT  of  the  target  depends  upon  the  emissivity,  target  temperature, 
target  size,  and  the  spatial  distribution  of  the  heat  (e.g.,  the  engine  area  will  be  the  hottest  area). 
The  equivalent  AT  is  calculated  so  that  the  MRT  can  be  used  to  describe  system  performance. 
Since  the  spatial  distribution  of  the  heat  and  the  temperature  can  change  from  test  to  test,  it  is 
easy  to  see  that  the  equivalent  AT  will  change.  Therefore,  it  becomes  exceedingly  difficult  to 
compare  field  data.  Conclusions  reached  about  the  effectiveness  of  smoke  may  result  from 
differences  in  the  detector  MRT  or  target  equivalent  AT  rather  than  the  optical  properties  of  the 
smoke. 


In  many  field  tests,  trained  observers  are  not  used.  The  relationship  between  a 
trained  observer  and  a  casual  observer  is  not  clear.  This  is  extremely  important  because  the 
systems  are  characterized  by  how  well  a  trained  observer  can  perceive  targets. 

VI.  TARGET  OBSCURATION. 

Consider,  now,  an  imaging  device  which  processes  the  detected  flux  into  various  gray 
levels  to  produce  an  image  on  a  TV  screen.  Let  us  assume  that  the  device  is  adjusted  so  that, 
before  the  smoke  is  introduced,  the  hot  target  will  appear  as  pure  white  and  the  background  will 
appear  as  black.  Let  us  assume  that  the  device  will  insert  1 0  gray  levels  between  these  points  so 
that  each  gray  level  is  given  by 


G  =  M 


!OT  ~ 
10 


(19) 


where  M  is  the  electronic  transfer  function  of  the  device. 

Before  any  further  analysis  is  possible,  the  exact  mode  of  operation  of  the  thermal 
imaging  system  must  be  known.  There  are  two  basic  modes.  In  the  first,  the  instrument  is 
adjusted  for  optimum  display  and  the  controls  are  not  further  adjusted,  even  after  the  smoke  is 
introduced.  In  the  second  mode,  the  internal  automatic  gain  control  (AGC)  automatically  adjusts 
the  gain  so  that  an  optimum  display  is  always  present.  The  first  mode  is  typical  of 
laboratory-type  systems  and  the  second  is  typical  of  military  imaging  (forward  looking 
infrared  —  FLIR)  systems.  Let  us  also  consider  the  case  where  the  smoke  is  homogeneous. 

Let  us  first  consider  a  device  which  will  not  be  readjusted.  There  are  three  possible 
conditions  which  will  obscure  the  target  (figure  16): 

CASE  A.  With  strong  attenuation  and  minimal  emission-reflection,  the  target 
signature  will  be  below  the  pure  black  level  and  the  entire  image  on  the  screen  will  be  black. 

!ST  <  *OB  <2°) 
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Figure  16.  Three  Methods  of  Obscuring  a  Target  Detected  with  a  Thermal  Imaging  System 

With  no  smoke,  the  target  is  pure  white  and  the  background  is  pure  black.  In 
A,  smoke  reduced  the  target  signature  below  the  black  level.  In  B,  the  smoke 
emission-reflection  raised  the  background  into  the  white  level.  In  C, 
emission-reflection  combined  with  attenuation  put  the  target-background 
intensity  differential  into  a  gray  level. 

which  results  in 


yWa' 


e-pR-p  +  La 


')dX</S(raeBRB  +  La)dX 


(21) 


If  La'  «  La,  then 


f$T  sTaeTRTdX  ^ f$T  aeBRBdX 


(22) 


Using  the  linearization  approximation  (equation  17)  and  letting  the  emissivities  be 
wavelength  independent  and  using  equation  10,  we  have 


<T>< 


eT  yl  +  a  AT  J 


(23) 


This  is  plotted  in  figure  17  for  both  the  3-  to  5-/im  and  8-  to  14-pm  regions.  The  graph  shows 
that,  with  minimal  emission-reflection,  the  transmission  can  be  quite  high  and  still  drive  the 
target  signature  below  the  background. 

For  the  system  that  is  not  readjusted,  the  reason  a  relatively  high  transmission  will 
obscure  a  target  is  that  the  flux  differential  between  the  target  and  background  is  very  small. 
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BACKGROUND— TARGET  TEMPERATURE  DIFFERENTIAL, 

Figure  17.  Average  Transmission  Necessary  to  Reduce  the  Target 
Signature  below  the  Black  Level  (Case  A) 


CASE  B.  With  strong  emission-reflection,  the  apparent  background  signal  (cloud 
plus  background)  will  increase  to  the  pure  white  level  and  the  entire  screen  will  be  white 


lSB  >  taT 


(24) 


which  results  in 


(25) 
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If  we  equate  the  emission-reflection  flux  to  some  equivalent  blackbody  temperature 
ATg  above  the  background  with  emissivity  eg  and  if  the  emissivities  are  wavelength  independent, 
then 


ATS>' 


( 1  +  a  AT)  e-p  -  <r>  eg  -  1 


esa 


(27) 


If  we  allow  ep  =  eg  =  1 ,  then 


AT, 


(28) 


For  es  =  1,  this  is  plotted  in  figure  18  for  both  the  3-  to  5-/im  and  8-  to  14-/im  systems.  The 
negative  ATg  indicates  that  the  cloud  can  be  colder  than  the  background  to  drive  the  background 
plus  cloud  signal  into  the  pure  white  range.  This  is  rather  interesting  for  it  indicates  that  it  does 
not  require  much  energy  (flux)  to  saturate  a  system  without  AGC. 
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Figure  18.  Required  Emission-Reflection  Temperature  Necessary  to 
Saturate  a  Thermal  Imaging  System  (Case  B) 
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CASE  C.  The  relationship  between  emission-reflection  and  attenuation  causes  the 
image  to  disappear  into  a  gray  level  and  then  the  entire  screen  will  be  gray.  Then 


!ST  ~  ]SB  < 


krr  -  (qb 

10 


(29) 


where 


!OB  <  ]SB 


]ST  ^  ]OT 

Comparing  this  to  equation  16,  we  see  that 


<r><To 


(30) 


This  is  an  unexpected  result!  It  indicates  that  the  observer  will  be  unable  to  see  his  target 
every  time  the  spectral  averaged  transmission  drops  below  10%  for  a  10  gray  level  system 
provided  that  he  does  not  readjust  his  thermal  imaging  system.  If  the  system  is  not  optimized, 
i.e.,  both  the  background  and  target  are  in  the  gray  levels,  then  the  relationship  presented  above 
is  relaxed  in  the  sense  that  a  higher  transmission  smoke  will  produce  the  same  effect.  Thus,  the 
three  cases  presented  are  “worst”  cases. 

Consider,  now,  the  device  which  has  an  AGC  or  one  that  is  readjusted  for  optimum 
display  after  the  smoke  is  present.  The  only  inequality  that  exists  is  when  the  attenuation  and 
emission-reflection  combined  produces  a  signal  differential  between  the  target  and  background 
that  is  below  the  equivalent  noise. 


]ST  '  JSB  <  lN 

If  we  represent  this  current  as  an  equivalent  blackbody  Rj^j  with  emissivity  e  =  1 ,  then 


(31) 


•B  ^B  ^^SRj^dX 


(32) 


The  equivalent  input  noise  is  defined  as  a  small  increment  about  the  background 
temperature  so  that 


RN  -  a  ATn  Rb 


(33) 
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Note  that  Rj^  is  defined  in  the  laboratory  and,  therefore,  the  atmospheric 
transmission  does  not  appear  in  the  right-hand  side  of  equation  33.  If  ej  =  eg  =  1,  then 


atn 

dX  =  <rs">  <  ~~~~ 
s  AT 


(34) 


which  is  the  well-known  equation  that  indicates  the  spectral  average  smoke  transmission  must  be 
sufficiently  small  to  reduce  the  target-background  differential  below  the  MRT  of  the  imaging 
system.  Note  also  that  this  average  <r">  is  different  from  that  defined  by  equation  1 9. 

For  both  systems  (e.g.,  with  AGC  and  without),  we  have  assumed  that  the 
emissivities  were  equal  to  one  in  several  places  and  that  the  blackbody  flux  was  wavelength 
independent  and  proportional  only  to  the  temperature.  These  assumptions  permitted  us  to  plot 
certain  parameters  and  to  draw  certain  conclusions.  In  practice,  however,  none  of  the 
assumptions  are  truly  valid,  and,  therefore,  we  must  consider  the  conclusions  as  guidelines  rather 
than  as  facts. 

Because  of  all  the  problems  associated  with  thermal  imaging  systems,  they  should 
not  be  used  to  quantify  the  effectiveness  of  smoke.  Rather  they  should  be  used  solely  for 
qualitative  purposes. 

VII.  THERMAL  CLUTTER. 

Thermal  clutter  is  defined  as  many  hot  sources  within  the  thermal  imaging  system’s 
FOV  that  will  confuse  the  observer  so  that  he  will  be  unable  to  locate  the  true  target.  The 
emissive  properties,  as  well  as  the  angular  subtense  required  of  each  individual  source,  are  not 
known  at  all.  Within  a  smoke,  natural  wind  velocity  and  diffusion  can  move  these  spots  so  that 
the  observer  will  see  only  randomly  moving  speckles  on  the  screen.  Some  of  the  hot  spots  should 
probably  be  at  least  the  size  of  the  target.  Pattern  recognition  studies  may  suggest  the  size, 
temperature,  and  number  of  local  hot  spots  required.  It  is  unknown  how  these  hot  spots  may  be 
generated. 


It  has  been  shown  in  previous  sections  that  a  uniform  emissive-reflective  cloud  only 
adds  a  dc  component  to  the  detected  flux.  The  electronics  are  capable  of  eliminating  this  shift 
and,  therefore,  the  target  will  be  discernible  within  the  limitations  already  mentioned.  It  was 
assumed  that  the  system  was  adjusted  so  that  the  target  was  pure  white  and  the  background  was 
pure  black.  In  some  systems,  the  device  locates  the  hottest  portion  of  the  scene  (presumably  the 
target)  and  the  coldest  portion  (presumably  the  background)  and  inserts  the  1 0  gray  levels  in 
between.  If  the  emission  is  sufficiently  high,  the  detector-electronics  will  saturate  causing  either  a 
portion  or  the  entire  display  to  bloom.  Consider,  now,  a  situation  where  a  portion  of  a  hot 
cloud  or  decoy  is  within  the  FOV.  If  the  detector  does  not  saturate  and  if  the  hot  area  is  only 
in  front  of  the  background  and  not  in  front  of  the  target,  then  the  pure  white  level  is  given  by 

yS(Ta'Rs+La'  +  TsTaeBRB)dX  (35) 
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and  the  pure  black  level  is  given  by 


S(raeBRB  +  La)dX 


(36) 


If  10  gray  levels  are  introduced  between  the  white  and  black  levels,  then  the  target  will  disappear 
into  a  gray  level  when 

/ Sra(eTRT  -  eBRfi)dX  <  ^/s(ra'Rs  +  La'  +  WBRB  -  raeBRfi  -  La)dX  (37) 


^  ^a  ~  ^a»  eT  eB  ^  >  ra  ~  Ta’  the  cloud  has  an  emissivity  es,  reflection  can  be  neglected,  and 
the  cloud  is  ATg  above  ambient,  then 


ATs>^[IOaAT  +1  "<r>'es]  (38) 

In  the  limit  as  the  smoke  becomes  extremely  dense,  <t>  -*■  0  and.  eg  -*■  1  then 

ATS  >  10  AT  (39) 

which  is  an  intuitive  result.  The  smoke  temperature  differential  above  the  background  must  be  at 
least  ten  times  the  target  differential  for  the  target  to  disappear  into  a  gray  level.  This  is  what 
can  happen  when  a  hot  cloud  starts  to  drift  into  the  FOV  of  a  thermal  imaging  system. 

VIII.  CONCLUSIONS. 

The  various  difficulties  encountered  in  measuring  transmission  of  aerosols  have  been 
discussed.  Specifically,  scattered  light,  whether  single  or  multiple,  can  enter  a  large  FOV  detector 
and  give  an  artificially  high  transmission.  The  size  distribution  of  an  aerosol  is  also  directly 
related  to  the  transmission.  Any  dynamic  process,  such  as  sedimentation  and  coagulation,  can 
alter  the  distribution  and  thereby  alter  the  extinction  measurements.  For  aerosol  measurements 
to  be  reproducible,  it  is  necessary  to  generate  exactly  the  same  size  distribution  and 
concentration  each  time.  A  flow  chart  indicating  possible  errors  in  measurement  for  reproducible 
aerosols  is  shown  in  Figure  1 9. 

To  intelligently  analyze  data  from  broadband  detectors,  the  spectral  components  of 
the  target,  smoke,  intervening  atmosphere,  and  the  detector  must  be  known.  To  obtain 
consistent  results,  in  addition  to  the  aerosol  reproducibility  requirements,  the  target  temperature 
and  intervening  atmosphere  must  always  be  the  same.  For  small  pathlengths,  the  atmospheric 
transmission  is  near  1 00%;  however,  for  long  pathlengths,  the  transmission  depends  upon  relative 
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LASER  TRANSMISSOMETER 


Figure  1 9.  Flow  Chart  of  Possible  Experimental  Errors 
Encountered  with  Laser  Transmissometers 


humidity,  temperature,  and  atmospheric  constituents,  all  of  which  can  vary  on  an  hourly  basis.  If 
two  different  detectors  are  used,  the  spectral  responses  must  be  identical.  The  flow  chart 
illustrating  these  problems  is  shown  in  figure  20. 


CONTINUE  TO 
LASER  TRANSMISSOMETER 
FLOW  CHART 


Figure  20.  Flow  Chart  of  Possible  Experimental  Errors 
Encountered  with  Broadband  Detectors 

Finally,  for  target  obscuration  with  a  thermal  imaging  system  with  an  AGC, 
sufficient  smoke  must  be  present  to  reduce  the  target-background  temperature  differential  below 
the  MRT  of  the  system.  The  MRT  depends  upon  the  angular  subtense  of  the  target.  The 
problems  inherent  to  broadband  detectors  and  transmissometers  also  apply.  The  flow  chart  is 
shown  in  figure  21.  Note  that  only  trained  observers  should  be  used  if  the  only  information 
available  is  the  display  screen. 
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CONTINUE  TO  BROADBAND 
DETECTOR  FLOW  CHART 


Figure  2 1 .  Flow  Chart  of  Possible  Experimental  Errors 
Encountered  with  Thermal  Imaging  Systems 
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APPENDIX  A 


EXAMINATION  OF  THE  CORRELATION  BETWEEN  LABORATORY 

AND  FIELD  SMOKE  EXTINCTION  DATA 


G.  C.  Holst 
M.  E.  Milham 

Chemical  Systems  Laboratory 
Research  Division 

Aberdeen  Proving  Ground,  Maryland 


I.  INTRODUCTION, 


When  lasers  or  narrow  bandwidth  detectors  are  used  to  obtain  extinction  coefficients 


of  various  smokes,  there  exists  reasonable  agreement  between  laboratory  and  field  measurements. 
However,  when  broadband  detectors  are  used,  the  measured  integrated  transmission  depends 
upon  the  spectral  transmission  of  the  smoke,  target  emissive  characteristics,  detector  spectral 
response,  and  the  atmosphere  spectral  transmission.  The  back  calculation  to  extinction 
coefficient  from  transmission  data  can  lead  to  any  arbitrary  result.  Thus,  with  knowledge  of  the 
relationship  among  the  various  parameters,  one  can  “prove”  that  the  extinction  coefficient  is 
almost  any  prechosen  value  by  judiciously  selecting  a  detector  with  the  appropriate  spectral 
response  and  a  target  of  the  appropriate  temperature. 


The  detected  signal  is  given  by 


(A-l) 


where 


T^  =  atmospheric  transmission 


R-p  =  energy  emitted  by  the  target 


Tg  =  smoke  transmission 


S  =  spectral  response  of  the  sensor 


and  the  integration  is  taken  over  the  wavelength  region  of  interest.  In  general,  all  parameters  are 
a  function  of  wavelength.  During  experiments,  the  usual  measurement  is  to  ratio  the  signal  with 
smoke  to  that  without  smoke  and  call  this  ratio  the  transmission.  This  is  represented  as 


If  the  smoke  has  no  spectral  absorption  (e.g.,  the  transmission  is  independent  of 
wavelength),  then  Tg  can  be  taken  out  of  the  integral  so  that  T=Tg  which  is  the  desired  result. 
With  narrow  band  detectors  or  lasers,  the  wavelength  region  is  so  small  that  T=Tg. 

Given  the  transmission,  it  is  tempting  to  infer  an  extinction  coefficient  a,  by  using 

Beer’s  law 


1  1 


(A-3) 


where 


C  =  concentration 
L  =  pathlength 

It  has  been  previously  shown1  that  the  effective  extinction  coefficient  (back 
calculated  from  transmission  and  CL  values)  is  a  function  of  the  target  temperature  and  sensor 
response.  Those  calculations  were  obtained  for  typical  3-  to  5-  and  8-  to  1 2-fJ.m  detectors. 

The  present  work  compares  field  data  obtained  at  Dugway  Proving  Ground  during 
Smoke  Week  I  with  theoretical  calculations.  These  calculations  include  the  actual  sensor  system 
and  target  temperature  used  at  Dugway  Proving  Ground.  Furthermore,  the  atmospheric 
transmission  using  LOWTRAN  4  has  been  included. 

II.  CALCULATIONS. 

The  spectral  extinction  coefficient  of  white  phosphorus  (WP)  obtained  in  the 
laboratory  is  shown  in  figure  A-l.^  To  illustrate  the  effect  of  spectral  extinction,  the  transmission 
of  WP  with  a  600-m  atmospheric  pathlength  is  plotted  in  figure  A-2  as  a  function  of  CL.  It  is 
readily  apparent  from  this  figure  that  since  broadband  detectors  integrate  under  the  curve, 
different  spectral  response  will  yield  different  total  transmissions.  Using  equation  A-2,  the 
transmission  of  several  commercially  available  detectors  has  been  calculated  (table  A-l). 

It  is  precisely,  the  fact  that  the  transmission  is  very  high  in  the  1 1-  to  13-/nm  region 
that  it  was  suggested  that  CO 2  isotope  lasers  should  be  considered. ^  The  expected  laser 
transmission  for  WP  is  given  in  table  A-2. 
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EXTINCTION  COEFFICIENT 


WHITE  PHOSPHORUS 


Figure  A-l.  Spectral  Extinction  Coefficient  of  WP 


WHITE  PHOSPHORUS 


Figure  A-2.  Spectral  Transmission  of  WP  for  Various  CL  Values 
(Atmospheric  pathlength,  600  m) 
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Table  A-l.  Calculated  Transmission  of  WP  Smoke  Using  Various  Detectors 


Detector 

CL  =  5 

CL  =  10 

% 

Infrared  associates  HgCdTe 

31 

14 

Barnes  immersed  thermistor 

26 

10 

Barnes  pyroelectric 

30 

13 

Barnes  HgCdTe 

31 

14 

Honeywell  HgCdTe 

25 

9 

Table  A-2.  Transmission  of  Various  CO2  Laser  Lines  in  WP 


CO2  laser  line 

CL  =  5 

CL  =  10 

pm 

% 

9.6 

7.4 

.55 

10.6 

19 

3.6 

11.2 

40 

16 

11.8 

52 

27 

Now  that  the  problems  associated  with  broadband  detectors  have  been  identified, 
consider  Dugway  Proving  Ground  test  trial  DP1-002-T23,  6  October  1977.4  The  munition  was  a 
155-mm  WP  round.  This  test  produced  CL  values  up  to  27.  The  atmospheric  conditions  of 
6  October  1977  were  inputted  to  the  LOWTRAN  code  to  obtain  the  atmospheric  transmission. 
Detector  response  curves  (bolometer)  and  filter  transmission  curves  (7.63  -  11.37  pm)  were 
obtained  from  Dugway  Proving  Ground.  The  source  was  an  1800°K  blackbody.  The  effective 
extinction  coefficient  was  back  calculated  using  equations  A-2  and  A-3  by  first  assuming  a  CL  value 
and  then  calculating  a.  A  plot  was  generated  for  all  values  of  CL  (figure  A-3).  The  actual  Dugway 
Proving  Ground  data  are  also  plotted  in  figure  A-3  where  the  transmission  from  the  “9.75-pm” 
filter  and  the  actual  CL  values  obtained  from  the  aerosol  photometers  were  used  to  calculate  an 
effective  extinction  coefficient. 


Appendix  A 


42 


Appendix  A 


43 


III.  DISCUSSION. 


The  fact  that  the  effective  extinction  decreases  as  CL  increases  is  a  result  of  the  fact 
that  a  broadband  detector  is  not  ideal  when  smokes  exhibit  spectral  character.  If  the  smoke  had 
constant  transmission  as  a  function  of  wavelength,  then  a  would  be  independent  of  CL.  As  a 
further  check  on  the  sensitivity  of  a  to  spectrally  absorbing  smokes,  calculations  were  performed 
on  lamp  black.  As  seen  in  figure  A-4,  lamp  black  has  a  very  small  change  in  a  as  a  function  of 
wavelength.  It  may  be  loosely  said  that  it  is  essentially  independent  of  wavelength  when 
compared  to  WP.  The  results  of  these  calculations  are  shown  in  figure  A-5 .  We  see  that  even  with 
this  slight  change  in  a  there  is  also  a  dependency  of  apparent  extinction  upon  CL. 

The  experimental  and  theoretical  curves  in  figure  A-3  are  quite  similar  in  shape  but 
differ  by  a  constant  factor.  The  fact  that  the  curves  are  similar  in  shape  supports  the 
methodology  used  to  calculate  the  extinction  coefficient.  There  are  several  possible  causes  for 
the  absolute  difference.  The  first  most  likely  is  that  the  laboratory  measurements  were  not 
calibrated  properly.  However,  recent  independent  laboratory  data5  obtained  with  a  fast  scanning 
Michelson  interferometer  agree  quite  well  with  Chemical  Systems  Laboratory  data.  The  second 
most  likely  cause  may  be  the  method  by  which  Dugway  Proving  Ground  obtains  extinction 
coefficients.  Their  method^  requires  knowledge  of  the  yield  factor;  whereas,  the  laboratory 
method  does  not. 

It  is  known  that,  when  phosphorus  is  oxidized,  it  does  not  immediately  hydrate  to 
orthophosphoric  acid. 2  Therefore,  the  theoretical  yield  factor^  which  applies  only  to 
orthophosphoric  acid  may  be  as  much  as  a  factor  of  two  higher  than  the  yield  factor  obtained 
when  phosphorus  is  initially  hydrated  into  various  metaphosphoric  and  tetrametaphosphoric 
acids.  Although  the  yield  factor  has  not  been  definitively  researched,  laboratory  yield  factors  are 
much  lower  than  the  theoretical  limits.  If  a  lower  yield  factor  were  used  by  Dugway  Proving 
Ground,  then  there  would  be  excellent  laboratory-field  agreement. 

If  the  detector’s  bandwidth  was  small  and  near  the  edge  of  transmission  (e.g.,  7  to 
8  pm),  then  the  atmosphere  would  significantly  affect  the  measured  transmission.  For  the 
detector-filter  combination  used  by  Dugway  Proving  Ground,  calculations  show  that  different 
environments  (hot  dry  versus  cold  damp)  do  not  appreciably  affect  the  measured  transmission. 

IV.  CONCLUSION  AND  RECOMMENDATION. 

» 

Broadband  detectors  are  not  ideal  when  measuring  transmission  of  smokes  unless  the 
detector  used  exactly  matches  a  fieldable  device  such  as  a  forward  looking  infrared  system. 
Because  of  the  dependence  of  the  extinction  coefficient  on  CL,  ratios  should  not  be  employed 
when  comparing  data  at  various  wavelengths. 

The  yield  factor  can  significantly  affect  the  Dugway  Proving  Ground  method  of 
obtaining  extinction  coefficients.  It  is  recommended  that  a  definitive  experiment  be  carried  out 
to  identify  each  intermediate  phosphoric  acid  species  and  to  determine  the  yield  factor.  Since 
the  species  evolve  in  time,  the  yield  factor  is  probably  a  function  of  time. 
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LAMPBLACK 


Figure  A-4.  Spectral  Extinction  Coefficient  of  Lampblack 
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Figure  A-5.  Calculated  Effective  Extinction  Coefficient  for  the 
“9.75-jim”  Detector  as  a  Function  of  CL 
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APPENDIX  B 


LINEARIZATION  OF  THE  BLACKBODY  CURVE 

As  was  mentioned  in  section  IV  and  shown  in  figure  14,  the  output  of  a  broadband 
imaging  system  is  calibrated  against  the  temperature  of  the  blackbody  source.  The  detector, 
however,  is  sensitive  to  the  flux  rather  than  the  temperature.  The  flux  is  given  by  Planck’s  law 

R(X,T)  =  KX‘5[EXP(B/XT)  -  1]_1  B-l 

where  K  =  constant,  B  =  14388  pm-°K,  X  =  wavelength,  and  T  =  the  temperature.  With  the  Wien 
approximation,  it  can  readily  be  shown  that  the  first  two  terms  of  the  Taylor  series  expansion 
are 


R(X,T)  *  R(X0,T0) 


B 

+ - ^  AT 

W 


(B-2) 


where  Tq  is  the  background  temperature  taken  as  Tq=  293°K,  Xq  is  a  wavelength  midpoint 
within  the  spectral  range  [Xq  =  4  pm  for  the  3-  to  5-pm  range,  and  Xq  =  1 1  pm  for  the  8-  to 
14-pm  range],  and  Tq  +  AT  is  the  temperature  of  the  source.  Therefore,  in  this  representation, 
AT  is  the  target-background  temperature  differential.  R(X,T)  is  independent  of  wavelength  and  is 
a  function  of  only  AT. 

The  error  associated  with  this  expansion  is  given  as 


R(X0,T0)  1  + 


E  =  1  - 


B 


XnT 


AT 


0*0 


R(X,T) 


(B-3) 


and  is  shown  in  figure  B-l.  The  errors  associated  with  the  3-  to  5-pm  region  are  extremely  large 
suggesting  that  this  representation  is  not  valid. 

The  reason  that  thermal  imaging  systems  can  be  linearized  is  that  the  output  is 
integrated  over  the  entire  spectral  range.  There  are  equal  positive  errors  and  negative  errors  in  the 
Taylor  series  expansion  and  these  tend  to  cancel.  To  prove  that  this  is  so,  the  intensity 

I  =  J  SRdX  (B-4) 

was  calculated  using  equation  B-l  (the  exact  solution)  and  with  equation  B-2  (the  approximate 
solution)  for  a  typical  InAs  (3-  to  5-pm)  detector  and  for  a  typical  HgCdTe  (8-  to  14-pm) 
detector.  The  relative  intensities  are  shown  in  figure  B-2.  Thus  the  integrated  output  is 
proportional  to  AT. 
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Note  the  different  ordinate  scales. 
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Figure  B-2.  Relative  Output  of  Two  Detectors  in  Two  Spectral  Ranges 

The  outputs  have  been  scaled  to  show  the 
linearity  between  output  and  AT. 


Appendix  B 


51 


DISTRIBUTION  LIST  6 


Names 

CHEMICAL  SYSTEMS  LABORATORY 

SAFETY  OFFICE 
Attn:  DRDAR-CLF 

AUTHOR’S  COPIES:  Research  Division 

USA  BIOMEDICAL  LABORATORY 
Attn:  SGRD-UV 

DEVELOPMENTAL  SUPPORT  DIVISION 
Attn:  DRDAR-CLJ-R 
Attn:  DRDAR-CLJ-L 
Attn:  DRDAR-CLJ-M 
Attn:  DRDAR-CLJ-P 

MUNITIONS  DIVISION 
Attn.  DRDAR-CLN 

RESEARCH  DIVISION 
Attn:  DRDAR-CLB 
Attn:  DRDAR-CLB-B 
Attn:  DRDAR-CLB-C 
Attn:  DRDAR-CLB-P 
Attn:  DRDAR-CLB-R 
Attn:  DRDAR-CLB-T 
Attn:  DRDAR-CLB-TE 

SYSTEMS  DEVELOPMENT  DIVISION 
Attn:  DRDAR-CLY-A 
Attn:  DRDAR-CLY-R 

DEPARTMENT  OF  DEFENSE 

Defense  Technical  Information  Center 
Attn:  DTIC-DDA-2 
Cameron  Station,  Bldg.  5 
Alexandria,  VA  22314 

Director 

Defense  Intelligence  Agency 
Attn:  DB-4G1 
Washington,  DC  20301 

DEPARTMENT  OF  THE  ARMY 

HQDA  (DAMO-SSC) 

WASH  DC  20310 

Director 

Defense  Civil  Preparedness  Agency 
Attn:  PO(DC) 

Washington,  DC  20301 

CINCUSAREUR 

Attn:  AEAGC-RSI 
APO  New  York  09403 

Deputy  Chief  of  Staff  for  Research, 
Development  &  Acquisition 
Attn:  DAMA-CSM-CM 
Attn:  DAMA-ARZ-D 
Washington,  DC  20310 


Copies 

I 

50 


12 


Names  Copies 

US  Army  Research  and  Standardization 

Group  (Europe)  1 

Attn:  Chief,  Chemistry  Branch 
Box  65,  FPO  New  York  09510 

HQDA  (DAMI-FIT)  I 

WASH,  DC  20310 

Commander 

HQ  US  Army  Medical  Command,  Europe 

Attn:  AEMPM  I 

APO  New  York  09403 

US  ARMY  HEALTH  SERVICE  COMMAND 

Superintendent 
Academy  of  Health  Sciences 
US  Army 

Attn:  HSA-CDC  1 

Attn:  HSA-IHE  I 

Fort  Sam  Houston,  TX  78234 

US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Attn:  DRCLDC  1 

Attn:  DRCSF-P  1 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 

Commander 

US  Army  Toxic  &  Hazardous  Materials  Agency 

Attn:  DRXTH-TS  2 

Aberdeen  Proving  Ground,  MD  21010 

Project  Manager  Smoke/Obscurance 

Attn:  DRCPM-SMK-M  1 

Aberdeen  Proving  Ground,  MD  21005 

Human  Engineering  Laboratory  HFE  Detachment 

Attn:  DRXHE-APG  1 

Building  520 

Aberdeen  Proving  Ground,  MD  21005 
Commander 

US  Army  Foreign  Science  &  Technology  Center 

Attn:  DRXST-MT3  1 

220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

US  Army  Materiel  Systems  Analysis  Activity 

Attn:  DRXSY-MP  1 

Attn:  DRXSY-T  (Mr.  Kaste)  1 

Aberdeen  Proving  Ground,  MD  21005 

Director 

DARCOM  Field  Safety  Activity 

Attn:  DRXOS-C  1 

Charlestown,  IN  47111 


53 


DISTRIBUTION  LIST  6  (Contd) 


Names 

Copies 

Names 

Copies 

Commander 

US  ARMY  TRAINING  &  DOCTRINE  COMMAND 

US  Army  Armament  Research  and  Development  Command 

Attn:  DRDAR-AC 

1 

Commandant 

Attn:  DRDAR-LCA 

1 

US  Army  Infantry  School 

Attn:  DRDAR-LCU 

1 

Attn:  NBC  Division 

1 

Attn:  DRDAR-SCA-C 

1 

Fort  Benning,  GA  31905 

Attn:  DRDAR-SCP-A 

1 

Attn:  DRDAR-SCW 

1 

Commandant 

Attn:  DRDAR-SER 

1 

US  Army  Missile  &  Munitions  Center  &  School 

Attn:  DRDAR-TSS 

2 

Attn:  ATSK-CD-MD 

1 

Dover,  NJ  07801 

Attn:  ATSK-DT-MU-EOD 

Redstone  Arsenal,  AL  35809 

1 

Commander 

Armament  Concepts  Office,  Weapon  Systems  Concepts  Team 

Commandant 

Attn:  DRDAR-ACW 

1 

US  Army  Military  Police  School/Training  Center 

Aberdeen  Proving  Ground,  MD  21010 

Attn:  ATZN-CDM 

1 

Attn:  ATZN-TDP-C 

4 

Director 

Ballistic  Research  Laboratory 

Fort  McClellan,  AL  36205 

Attn:  DRDAR-TSB-S 

2 

Commander 

Building  305 

US  Army  Infantry  Center 

Aberdeen  Proving  Ground,  MD  21005 

Attn:  ATSH-CD-MS-C 

1 

Fort  Benning,  GA  31905 

CDR,  APG 

USA  ARRADCOM 

Commandant 

Attn:  DRDAR-GCL 

1 

US  Army  Ordnance  &  Chemical  Center  &  School 

Aberdeen  Proving  Ground,  MD  21010 

Attn:  ATSL-CL-CD 

1 

Aberdeen  Proving  Ground,  MD  21005 

US  ARMY  ARMAMENT  MATERIEL  READINESS  COMMAND 

Commander 

Commander 

USA  Training  &  Doctrine  Command 

US  Army  Armament  Materiel  Readiness  Command 

Attn:  ATCD-Z 

1 

Attn'  DRSAR-LE 

1 

Fort  Monroe,  VA  23651 

Attn:  DRSAR-PDM 

1 

Attn:  DRSAR-SA 

1 

Commander 

Attn.  DRSAR-SF 

1 

USA  Combined  Arms  Combat  Development  Activity 

Rock  Island,  IL  61299 

Attn:  ATZLCA-CO 

Fort  Leavenworth,  KS  66027 

1 

CDR, APG 

USA  ARRCOM 

US  ARMY  TEST  &  EVALUATION  COMMAND 

Attn:  SARTE 

1 

Aberdeen  Proving  Ground,  MD  21010 

Commander 

US  Army  Test  &  Evaluation  Command 

Commander 

Attn:  DRSTE-FA 

1 

US  Army  Dugway  Proving  Ground 

Aberdeen  Proving  Ground,  MD  21005 

Attn:  Technical  Library,  Docu  Sect 

Dugway,  UT  84022 

1 

DEPARTMENT  OF  THE  NAVY 

Commander 

Commander 

Rocky  Mountain  Arsenal 

Naval  Explosive  Ordnance  Disposal  Facility 

Attn:  SARRM-QA 

1 

Attn:  Army  Chemical  Officer,  Code  604 

1 

Commerce  City,  CO  80022 

Indian  Head,  MD  20640 

Commander 

Commander 

Pine  Bluff  Arsenal 

Nuclear  Weapons  Training  Group,  Atlantic 

Attn:  SARPB-ETA 

1 

Naval  Air  Station 

Pine  Bluff,  AR  71611 

Attn:  Code  21 

1 

Norfolk,  VA  23511 

54 


DISTRIBUTION  LIST  6  (Contd) 


Names 


Copies 


Names 


Copies 


Chief,  Bureau  of  Medicine  &  Surgery 
Department  of  the  Navy 
Washington,  DC  20372 


ADDITIONAL  ADDRESSEES 


Commander 

USEUCOM 


Commander 

Naval  Weapons  Center 


Attn:  ECJ-O/LTC  James  H.  Alley 
APO,  NY  09128 


Attn:  A.  B.  Galloway /Code  3171 
China  Lake,  CA  93555 


Commanding  Officer 

Naval  Weapons  Support  Center 

Attn:  Code  5042 /Dr.  B.  E.  Douda  1 

Crane,  IN  47522 

US  MARINE  CORPS 

Director,  Development  Center 

Marine  Corps  Development  &  Education  Command 

Attn:  Fire  Power  Division  1 

Quantico,  VA  22134 

DEPARTMENT  OF  THE  AIR  FORCE 

HQ  Foreign  Technology  Division  (AFSC) 

Attn:  PDRR  1 

Wright-Patterson  AFB,  OH  45433 

Commander 

Aeronautical  Systems  Division 

Attn:  ASD/AELD  1 

Wright-Patterson  AFB,  OH  45433 

HQ,  USAF/SGES  *  1 

Bolling  AFB,  DC  20332 

HQ  AF1SC/SEV  1 

Norton  AFB,  CA  92409 

OUTSIDE  AGENCIES 

Battelle,  Columbus  Laboratories 

Attn:  TACTEC  1 

505  King  Avenue 
Columbus,  OH  43201 

Director  of  Toxicology  1 

National  Research  Council 
2101  Constitution  Ave,  NW 
Washington,  DC  20418 


55 


